Wigner function for twisted photons 
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A comprehensive theory of the Weyl-Wigner formalism for the canonical pair angle-angular momentum is 
presented, with special emphasis in the implications of rotational periodicity and angular-momentum discrete- 
ness. 
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I. INTRODUCTION 

A quantum system has a dynamical symmetry group G if 
its Hamiltonian is a function of the generators of G. In this 
case, the Hilbert space of the system spUts into a direct sum 
invariant subspaces (carriers of the irreducible representations 
of G) and the discussion of any physical property can be re- 
stricted to one of these subspaces yj). 

The existence of such a symmetry also allows for the ex- 
plicit construction of a phase space for the system as the 
coadjoint orbit associated with an irreducible representation 
of G iQIltl (in fact, it turns out to be a symplectic manifold). 
In consequence, to every operator on Hilbert space we can 
associate a function on phase space, opening the way to for- 
mally representing quantum mechanics as a statistical theory 
on classical phase space. Various aspects of this formalism for 
basic quantum systems have been developed by a number of 
authors 

There are, however, important differences with respect to a 
classical description. They come from the noncommuting na- 
ture of conjugate quantities, which precludes their simultane- 
ous precise measurement and, therefore, imposes a fundamen- 
tal limit to the accuracy with which we can determine a point 
in phase space. As a distinctive consequence of this, there 
is no unique rule by which we can associate a classical phase- 
space variable to a quantum operator and depending on the op- 
erator ordering, various functions can be defined. For exam- 
ple, the quantum state (i.e., the density matrix) of the system 
can be mapped onto a whole family of functions parametrized 
by a number s; the values +1, 0, and —1 corresponding to 
the Husimi Q, the Wigner W, and the Glauber-Sudarshan 
P functions, respectively. These phase-space functions are 
known as quasiprobability distributions, as in quantum me- 
chanics they play a role similar to that of genuine probability 
distributions in classical statistical mechanics (for reviews, see 
Refs. 116i,il7),a8.,.19J). 

Apart from the description of the harmonic oscillator (for 
which G is the Heisenberg-Weyl group and the correspond- 
ing phase space is the plane M^), this formalism has also 
been successfully applied to spin-like systems (or qubits in 
the modern parlance of quantum information), for which G is 
the group SU(2) and the phase space is the two-dimensional 
Bloch sphere. However, one can rightly argue that this Wigner 
function, although describing a discrete system, is not defined 



in a discrete phase space. In fact, the growing interest in quan- 
tum information has fueled the search for discrete phase-space 
counterparts of the Wigner function (see Ref. 120[| for a com- 
plete and up-to-date review). The main advantage of such a 
representation consists in that even states from different ir- 
reducible representations can be pictured on the same phase 
space, which is basically a direct product of two-dimensional 
discrete tori. 

There is still another "mixed" canonical pair: angle and an- 
gular momentum. Now, the symmetry group G is noncom- 
pact and can be taken as the two-dimensional Euclidean group 
E(2), whereas the associated phase space is the discrete cylin- 
der Z X 5i {Si denotes here the unit circle), since one of the 
variables is continuous and the other is discrete. Several inter- 
esting properties of a number of systems, such as molecular 
rotations, electron wave packets. Hall fluids, and light fields, 
to cite only a few examples, can be described in terms of this 
symmetry group 12111 . In quantum optics, it is the basic tool 
to deal with the orbital angular momentum of the so-called 
twisted photons ll22[|23ll . which have been proposed for appli- 
cations in quantum experiments ll24ll . 

The construction of a proper Wigner function for this case 
is still under discussion. Although some interesting attempts 
have been published lEsl l26l EtIi . they seem of difficult 
application to practical problems. Quite interesting group- 
theoretical approaches to this problem can be also found in 
Refs. 1281 l29ll . In this paper, we approach this interesting 
problem from the perspective of finite-dimensional systems 
and construct a bona fide Wigner function that fulfills all the 
reasonable requirements and is easy to handle and to inter- 
pret. We also discuss its applications to some relevant quan- 
tum states. 



II. WIGNER FUNCTION FOR POSITION-MOMENTUM 

In this section we briefly recall the relevant structures 
needed to set up the Wigner function for Cartesian quantum 
mechanics. This is to facilitate comparison with the angular 
case later on. For simplicity, we choose one degree of freedom 
only, so the associated phase space is the plane R^. 

The canonical Heisenberg commutation relations between 
Hermitian coordinate and momentum operators q and p are (in 
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units H = 1) 

[q,p]=t, (2.1) 

SO that they are the generators of the Heisenberg-Weyl alge- 
bra. In the unitary Weyl form this is expressed as 

U{p)V{q)=V{q)U{p)e"'P, (2.2) 

where 

V{q) = cxp(-iqp) , U{p) ^ cxp(ipq) , (2.3) 

are the generators of translations in position and momentum, 
respectively. In the Cartesian case, these exponentials can be 
entangled to define a displacement operator 

b[q,p) = U[p)V{q)e-''^P'^ = exp[z(pg - qp)] , (2.4) 

with the parameters {q^p) labelling phase-space points. How- 
ever, this cannot be done for other canonical pairs, as we shall 
see. 

The displacement operators form a complete trace- 
orthonormal set (in the continuum sense) in the space of op- 
erators acting on Ti (the Hilbert space of square integrable 
functions on M): 

Tl[b{q,p)b\q',p')] = 2Tl5{q-q')5{p-p'). (2.5) 

NotethatZ)t(,7,p) = -p), while Z)(0, 0) = 1. 

The mapping of the density matrix g into a Wigner function 
defined on is established in a canonical way: 

W{q,p) = Tr[^ui(g,p)] , 

(2.6) 

Q=j7r-Tn w{q,p)W{q,p)dqdp, 
(27r)^ 

where the (Hermitian) Wigner kernel w (a particular instance 
of a Stratonovitch-Weyl quantizer) is the double Fourier trans- 
form of the displacement operator: 

'^(Q^p) ^ 777-^ [ cxp[-i{pq' ~ qp')]b{q',p')dq'dp' . 

(2.7) 

One can immediately check that the Wigner kernels are also a 
complete trace-orthonormal set. Furthermore, they transform 
properly under displacements 

w{q,p) - b{q,p) w{Q, 0) b\q,p) , (2.8) 

where 

w{Q,Q)^ I b{q,p)dqdp^2P, (2.9) 

and P is the parity operator 

The Wigner function in ( 12.6b fulfills all the basic properties 
required for any good probabilistic description. First, due to 
the Hermiticity of w{q,p), it is real for Hermitian operators. 



Second, on integrating W{q,p) over one variable, the proba- 
bility distribution of the conjugate variable is reproduced 

/ W{q,p)dp= {q\g\q) , / W{q,p) dq = {p\q\p) . 
Jm JR 

(2.10) 

Third, W{q,p) is covariant, which means that for the dis- 
placed state g' = b{qo,po) g&{qo,po), one has 

W^,{q,p) = Wi,iq-qo,p-po), (2.11) 

so that the Wigner function follows displacements rigidly 
without changing its form, reflecting the fact that physics 
should not depend on a certain choice of the origin. 

Finally, the overlap of two density operators is proportional 
to the integral of the associated Wigner functions: 

Tr(eie2)«/ Wi{q,p)W2{q,p)dqdp. (2.12) 

This property (often known as traciality) offers practical ad- 
vantages, since it allows one to predict the statistics of any 
outcome, once the Wigner function of the measured state is 
known. 



III. WIGNER FUNCTION FOR DISCRETE SYSTEMS 

Many quantum systems can be appropriately described in 
a finite-dimensional Hilbert space. The previous standard ap- 
proach can be extended to these discrete systems, since they 
do have a dynamical symmetry group. However, in a continu- 
ous Wigner function for these systems, there is a lot of infor- 
mation redundancy. The goal of this section is to carry out a 
non-redundant discrete phase-space analysis for this case. 

Let us consider a system living in a Hilbert space Hd, of 
dimension d (a qudit). It is useful to choose a computational 
basis |n) (n = 0, . . . ,d — 1) in Hd and introduce the basic 
operators OOll 

X\n) = \n + l), Z\n) = uj{n)\n) , (3.1) 

where addition and multiplication must be understood modulo 
d and, for simplicity, we use the notation 

w(m) = uj'"- = exp(i27rm/d) , (3.2) 

u! = exp{i2TT / d) being a dth root of the unity. The operators 
X and Z generate a group under multiplication known as the 
generahzed Pauli group Olll and obey 

ZX^LuXZ, (3.3) 

which is the finite-dimensional version of the Weyl form ( 12.21 ) 
of the commutation relations. 

The monomials {Z*'X^} {k,l — . . . , d — 1) form a ba- 
sis in the space of all the operators acting in Tid Il32ll . It seems 
then natural to introduce the unitary displacement operators 

b{k,l) = e"^^^^^^Z^X\ (3.4) 
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where 4>{k, I) is a phase. The unitarity condition imposes that 

(j){k,l) + <p{-k,-l)^-'^kl. (3.5) 

Different choices have been analyzed in the Hterature ll33ll : 
one of special relevance is 

(j){k,l) = ^2-Hl, (3.6) 

where 2^^ is the multiplicative inverse of 2 in when d is 
prime and 2^^ = 1/2 for nonprime dimensions. 

In this way, we have got a discrete phase space of the system 
as a d X d grid of points, in a such a way that the coordinate 
of each point (fc, I) define powers of Z ("position") and X 
("momentum") and the whole phase space is isomorphic to a 
discrete two-dimensional torus. 

The following mapping from the Hilbert space into the dis- 
crete phase space [equivalent to ( 12.61 )1 

W{k,l) = Ti[gw{k,l)] , 

(3.7) 

k,l 

is established in terms of the following (Hermitian) Wigner 
kernel 

w{k,l) = -^'^^Lij{kn ~ Im) D{m,n) , (3.8) 

which is normalized, satisfies the overlap condition 

Tr[u;(fc, l)wik', I')] = d5k,k' ki' . (3-9) 
and it is explicitly covariant: 

wik,l) = Dik,l)wiO,0) D^kJ) , (3.10) 

where 

u'(0,0) = ^^l)(fc,0- (3.11) 

It is interesting to note that the phase (13.6b for prime dimen- 
sions leads to w{Q, 0) = P, P being the parity operator In 
view of these properties, one can easily conclude that the cor- 
responding Wigner function W{k, I) fulfills properties fully 
analogous as those for the continuous harmonic oscillator. 

IV. WIGNER FUNCTION FOR ANGLE-ANGULAR 
MOMENTUM 

In this section, we consider the conjugate pair angle and 
angular momentum. To avoid the difficulties linked with pe- 
riodicity, the simplest solution ll34l [35i [36ll is to adopt two 
angular coordinates, such as, e.g., cosine and sine, we shall 
denote by C and S to make no further assumptions about the 



angle itself. One can concisely condense all this information 
using the complex exponential of the angle E = C + iS, 
which satisfies the commutation relation 

[E,L] = E, (4.1) 

or, equivalently, 

(4.2) 

[C,S] = 0. 

In mathematical terms, this defines the Lie algebra of the two- 
dimensional Euclidean group E(2). Note also, that from the 
Baker-Campbell-Hausdorff formula, one gets 

e-"f'^E = e'"^ Ee~''l'^ , (4.3) 

which is the unitary Weyl form of ( 14.11 ). 

The action of E on the angular momentum basis is 

E\£)^\£-l), (4.4) 

and, since the integer £ runs from — oo to +oo, E is a unitary 
operator whose normalized eigenvectors 

form a complete basis 

(010') = ^ 5(0 - 0' - 2£7r) = 52.(0 - 0') , (4.6) 

where 62-^ represents the periodic delta function (or Dirac 
comb) of period 2tt. 

As anticipated in the Introduction, the phase space is now 
the semi-discrete cylinder x Si. Following the ideas of 
SecHni a displacement operator can be introduced as 

D{£^ 0) = e^"^^''^^ E-'^e-"^^ , (4.7) 

where a{£, 0) is a phase to be specified. Note that here there is 
no possibility to rewrite Eq. ( 14. 7t as an entangled exponential, 
since the action of the operator to be exponentiated would not 
be well defined. The requirement of unitarity imposes now 

a{£,(j)) + a{-£,-(j)) = £(/}. (4.8) 

As desired, the displacement operators form a complete 
trace-orthonormal set; 

Tr[7^(£,0)L)t(£'^^')] ^27r5«, 52.(0-0'), (4.9) 

whose resemblance with relation ( 12.5b is evident. 
We can introduce then the canonical mapping 

W{£,<l>) = Tr[gw{£,^)], 

(4.10) 
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where the Wigner kernel w is defined, in close analogy to the 
previous cases, as 
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4{e'<j) - e(j)')]D{£', (j)') d4>' . 



(4.11) 

The set of Wigner kernels constitutes a complete orthogo- 
nal Hermitian operator basis. In addition, they are explicitly 
covariant: 



w 



(€,0) = £)(A0)u;(0, 0)-Dt (4.12) 



with 



^i(0,0) = -i-2 5^ / mA)dc^, (4.13) 

although the interpretation of t&(0, 0) as the parity on the 
cylinder is problematic. 

All these properties automatically guarantee that we have 
indeed a well-behaved Wigner function for this canonical pair. 



V. EXAMPLES 

To work out the explicit form of the Wigner function for 
a given state, one first needs to specify the phase a{£,4') in 
Eq. ( I4.8l l. For convenience, in this paper the choice 



(5.1) 



shall be used, as it is linear in both arguments, and it appears to 
be the simplest function fulfilhng the unitarity condition and 
the periodicity in 4> lf2lll . 

In this case, the Wigner kernel (14.11b becomes 



- y 



After some manipulations, we obtain 



(5.2) 



— y 



( 1)^ ^-(2<;'+i)i0 



27r2 ^ I" + 1/2 

X + r + + £"-£'1, (5.3) 

which coincides with the kernel derived by Plebanski and 
coworkers |l29'] in the context of deformation quantization. 

Note that ( 15.31 ) splits into "even" and "odd" parts, depend- 
ing on whether the matrix elements qui ~ (t\Q\t') have ^ ± £' 
even (first sum) or odd (second sum). 

For an angular momentum eigenstate one immediately 
gets 



(5.4) 



which is quite reasonable in this case: it is flat in (\) and the 
integral over the whole phase space gives the unity, reflecting 
the normalization of |£o). 



For an angle eigenstate |0o), one has 



(5.5) 



Now, the Wigner function is flat in the conjugate variable ^, 
and thus, the integral over the whole phase space diverges, 
which is a consequence of the fact that the state |(/)o) is not 
normalizable. 

The coherent states (pp) (parametrized by points on the 
cylinder) introduced in Ref. IstIi (see also Refs. ifJsl [39ll for 
a detailed discussion of the properties of these relevant states) 
are characterized by 
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0, (Po) 



(014, (Ao) 



Jf.oi'p-'t'a) 



g-if0o ^-(l-loY/2 



(5.6) 



where -dz denotes the third Jacobi theta function 14C 
The Wigner function for the state \Iq,4>q) splits as 



(+) 



The "even" part turns out to be 



(-) / 



(5.7) 



(+) 



27r^3 (0|i) 



(5.8) 

This seems a sensible result, since it is a discrete Gaussian in 
the variable £, and for the continuous angle it is a Jacobi 
theta function, which plays the role of the Gaussian for circu- 
lar statistics |)4Tll. However, the "odd" contribution spoils this 
simple picture: 



(-) ( 

\f-aAoV 



gj(0-0o)-l/2 



?' + 4-^+l/2 



(5.9) 



In Fig. [T] the Wigner function for the coherent state |£o = 
0, 4>o = 0) is plotted on the discrete cylinder. A pronounced 
peak at (/) = for I — Q and slightly smaller ones for ^ = ±1 
can be observed. The associated marginal distributions [ob- 
tained from Eq. ( 15. 7t by integrating over or by summing 
over respectively] are also plotted. They are strictly posi- 
tive, as correspond to true probability distributions. 

For quantitative comparisons, however, sometimes it may 
be convenient to "cut" this cylindrical plot along a line 
(/)=constant and unwrap it. This is shown in Fig. |2] Here, the 
range of I is from -4 to 4, while the angle is plotted between 

— TT to TT. 



FIG. 1 : Plot of the Wigner function for a coherent state with £o — 
and (po = 0. The cylinder extends vertically from ^ = — 4 to ^ = +4. 
The two corresponding marginal distributions are shown. 

A closer look at these figures reveals also a remarkable fact: 
for values close to (j) ~ ±tt and £ — ±1, the Wigner function 
takes negative values. Actually, a numeric analysis suggests 
the existence of negativities close to = ±7r for any odd 
value of £. 

As our last example, we address the superposition 

l*) = ^(l^i)+e'^14)) (5.10) 
of two angular-momentum eigenstates with a relative phase 



FIG. 3: Plot and marginal distributions of the Wigner function for an 
even superposition \£i +g £2) with £1,2 = ±3 for £ — —4 to £ = +4. 

gitpo analysis can be carried out for the superposition of 
any number of eigenstates, but dS.lOt is enough to display the 
relevant features. 

The Wigner function splits again; now the "even" part reads 

as 

+ 2Si,+e,,2i cos[0o + (^2 -^i)<^]}(5.11) 




FIG. 2: Unwrapped plot of the Wigner function for a coherent state 
with ^0 = and (j)o = 0. The plane extends from £ = —4 to £ = +4 
and from <p = —n to (f) = n. 



FIG. 4: Plot and marginal distributions of the Wigner function for an 
even superposition \£i +g £2) with £1 = 4, £2 = —3 for £ = —4 to 
£ = +5. 
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For the "odd" part, the diagonal contributions vanish, and one 
has 

Wl-^{£,4:) = ^COs[0o + (^2-^l)0] 
(_l)^+(^i+'^2-l)/2 

^ ^,+t,-2i ^^i+^.=odd,(5.12) 

where 6ti+i2=odd indicates that the sum is nonzero only when 
li + I2 is odd. 

In consequence, when |£i — €2 1 is odd, the interference term 
contains contributions for any £, damped as When |£i — 
^2! is an even number, the contribution ( 15.12b vanishes and 
we have three contributions: two flat slices coming from the 
states l^i) and 1^2) and an interference term located at ^ = 

(^l+^2)/2. 

These features are illustrated in Figs.[3]and|4l The state \^) 
is plotted for £2 = ~3 and ti ~ 'd and (Fig. [3]) and I2 = ~3 
and ^1=4 (Fig.|4]i. Changing the relative phase 0o results in 
a global rotation of the cylinder In can be observed in Fig. |4] 
that the two rings at £ = —3 and £ = 4 (as opposed to the 
rings at £ = ±3 in Fig.[3]l, are not flat in 0, but show a weak 
dependence on the angle due to the odd contributions added 



to the flat Kronecker deltas. 



VI. CONCLUDING REMARKS 

In summary, we have carried out a full program for a com- 
plete phase-space description in terms a Wigner function for 
the canonical pair angle-angular momentum. An experimen- 
tal demonstration in terms of optical beams is presently un- 
derway in our laboratory. 
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